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We calculate the power spectra of primordial curvature and isocurvature perturbations from 
a general two field inflation model at next-to-leading order correction in a slow-roll expansion. In 
particular we calculate the spectral indices to second order in slow-roll parameters. We show that the 
cross-correlation of the curvature and isocurvature perturbations at the time of Hubble-exit during 
inflation is non-zero at first-order in slow-roll parameters. We apply our results to different classes 
of inflation, including inflaton and curvaton scenarios. The spectrum of primordial gravitational 
■ waves, curvature and isocurvature perturbations obey generalised consistency relations in two-field 

inflation models. We give the first two consistency relations in an infinite hierarchy. 

o ' 

(N) ■ PACS numbers: 98.80.Cq astro-ph/0605679 



in 



i 

o 
U 



I. INTRODUCTION 



During inflation the vacuum fluctuations of light scalar fields grow into super-Hubblc density perturbations which 
are believed to be the origin of the structure seen in the universe today. Single-field slow roll inflation predicts a 
C*~) ■ nearly scale invariant spectrum of gaussian curvature perturbations. However from a particle physics point of view it 
is natural to expect there to be more then one scalar field rolling during inflation. During multiple field inflation a 

■ spectrum of isocurvature as well as curvature perturbations can be generated and the two may be correlated [1] . 
In this paper we consider the case of two field inflation with an arbitrary potential and arbitrary background tra- 
jectory. Following Rcf. [2, 3] we introduce arbitrary transfer functions to parameterise the generation of primordial 
curvature and isocurvature perturbations from isocurvature perturbations during inflation. We present the power 

\Q , spectra of the isocurvature and curvature perturbations to first-order in slow roll as well as the cross-correlation, 

■ presenting the results in terms of the slow-roll parameters and the (in principle) observable correlation angle. We are 
*"fH then able to calculate the tilts of all three power spectra to second-order in slow-roll parameters (see also Ref. [4]). 

Q. This generalises the result of Stewart and Lyth [5] for adiabatic perturbations from single-field inflation. For an alter- 
native approach to calculating the primordial curvature perturbation in multi-field models (but not the isocurvature 
perturbation), based on the SN- formalism [6], see Ref. [7, 8]. 

In section II we calculate the scalar field perturbations during inflation, first performing a global rotation in field 
space to a basis where the perturbations are uncorrelated at Hubble exit. We then rotate to the local (instantaneous) 
basis of adiabatic and isocurvature perturbations [4, 9, 10] and show that they are correlated at first-order in slow roll 
at Hubble exit. The only exception is the case of a straight background trajectory in field space, e.g. for a symmetric 
potential [11], when the local rotation is in fact a constant (global) rotation which coincides with the uncorrelated 
basis and hence the curvature and isocurvature perturbations are uncorrelated at Hubble exit. 

In section III we change variables to the dimensionless curvature and isocurvature perturbations which are then 
related to the primordial curvature and isocurvature perturbations after inflation, and introduce the dimensionless 
correlation angle. The curvature and isocurvature perturbations are evolved on supcrhorizon scales by introducing two 
arbitrary transfer functions, which parameterise our ignorance of the detailed physics after Hubble exit and through 
the end of inflation and reheating. We present results for the power spectra and tilts including the next to leading 
order terms and also the running at leading order, which is second-order in slow roll. We also briefly present the 
power spectrum, tilt and running for tensor perturbations, which are unchanged from the single-field case [12]. In 
section IV we simplify the results to the interesting special cases of a straight background trajectory, which includes 
assisted inflation with exponential potentials [13-15], the curvaton scenario [16], and the inflaton scenario allowing 
for an additional uncorrelated isocurvature field present during inflation. We conclude in section V. 

Appendices contain definitions of all the slow-roll parameters used and the relations between these parameters when 
defined in terms of the Hubble parameter or the potential, their derivatives and further details of the calculation of 
the tilts. Throughout this paper we use the notation that results accurate at first order in slow roll are denoted by 
~, while equality at second order in slow roll is denoted by =. 
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II. INITIAL POWER SPECTRA 



A. Background equations of motion 



We consider two scalar fields with Lagrange density 

i 2 



(i) 



7=1 



We thus allow for an arbitrary interaction potential, V(</>i,(fo), but for simplicity consider only canonical kinetic 
terms. It should be possible to generalise our results to non- minimal kinetic terms, see for example [4, 10]. 
The background equations of motion are 



b! + 3i707 + Vi = 0, 



H 



8ttG 



v + \v 



where Vj = dV/d<pi. To solve these equations we use the slow roll approximation to rewrite them as 



8ttG, 



1 + 



ii 



(2) 
(3) 



(4) 



where 5f and e H are slow-roll parameters as defined in appendix A, which we assume to be small. We make the 
standard slow roll assumption that the time derivative of the slow-roll parameters are higher order in slow roll. 



B. The perturbation equations 



The equations of motion for the perturbed fields in the spatially flat gauge are [17] 



5<j>j + ZHSfa + —5<t>i + Y. 



Vu- 



8irG d 



a 3 dt\H 



'JI9.1 



Hi 



(5) 



where k is the comoving wavenumber. To simplify the above equations we change variables to uj = aS4>i and to 
conformal time r, defined by dr = dt/a. To lowest order in slow roll, r ~ —(1 + e)/(aH), the equations of motion 
simplifies to 



fc2 - Z2 ) U I = Z2 J2 MljUj 



(G) 



where the prime denotes a derivative with respect to conformal time and the interaction matrix is 



Vl2 



2ef 2 - Vi2 e H + 2e£ - ry 22 



(7) 



Note that the superscript H denotes a slow-roll parameter defined in terms of the Hubble parameter, while all other 
slow- roll parameters are defined in terms of the potential (see appendix A). 

The two differential equations for ui are coupled but they can be decoupled at Hubble-exit by performing a rotation 
to diagonalise the interaction matrix (7) to first order in slow roll. For previous discussions on how to solve the coupled 
equations (6) see, for example, [10, 18]. The rotation matrix is given by 



U = 



where the rotation angle is defined by 

tan 26 = 2 



cos 9 — sin 
sin cos 



2e? 2 - JJia 



2(efi - eg) - (7?ii - r? 22 ) 



(8) 



(9) 
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The subscript * refers to a quantity evaluated at horizon crossing, k = a*H*. For a given wavenumber k, <d is a 
constant, but there is a dependence since the slow roll parameters are evaluated at horizon crossing and each wave 
mode has a different horizon crossing time. Note that since the rotation angle 8 is defined in terms of a ratio of slow- 
roll parameters it is in general not small. If all of the slow-roll parameters become arbitrarily small Eq. (9) becomes 
ambiguous, because the equations of motion (6) are then decoupled to first order in any basis and any rotation will 
work. The rotation matrix (8) satisfies 

C/ t Aft/ = diag(Ai,A 2 ), (10) 

where 



(11) 



Ai,2 *\ 4e H - fan + r? 22 ) ± yf[2(e» - eg) - (rj u - + 4(2ef 2 - Vl2 f 

Denoting the decoupled variables by uj, we have 

ui = ]T U ijV j . (12) 

Therefore left-multiplying (6) by and rewriting in the vi basis we have 

v" I + (k> l <-> 1 

where 



* -72(^-4) K-O, (13) 



Mi = \ + A/ . (14) 

This is the same form for the equation of motion as in the single field case [5] except that the /i/ depend on slow-roll 
parameters relating to both the original fields (f>\ and (f>2 ■ 

Working at first order in slow roll, so that we can treat the slow-roll parameters and hence \xj as constant, the 
solution is 

Vl = ^e^^+ 1 /2W2 ( _ T) i/2 jff a) ( _ fcr)e/(fc) (15) 

where H^~} '(x) is a Hankel function of the first kind of order \xi and the ei(k) are independent unit Gaussian random 
variables [9] satisfying 

(e/(k)) = , (ej(k)e* (kO) = S L j5 3 (k - k') . (16) 
The constant in front of the Hankel function solution for m comes from imposing the early-time boundary condition 

vj — ► -^=e~ ikT eAk) as - kr — > 00 . (17) 
V2fc 



The corresponding late-time behavior (fcr — > 0) is 

T(3/2)V2fc 7 v ; 



ie iA lW /2^ + c\i)-^={-kr)- 1 - Xt ei(k) (18) 
v2fc 



where C = 2 — log 2 — 7 w 0.7296 and 7 ~ 0.5772 is the Euler-Mascheroni constant. Note however that although the 
change in first-order slow-roll parameters and hence \ii in one Hubble time is second order in slow roll, see appendix 
B, the time variation of the slow-roll parameters cannot be neglected over many Hubble times. Thus we can only 
reliably use the late-time solution (18) for a few Hubble times after Hubble-exit (k = a*H*). We need to perform 
another rotation in field space to accurately track the evolution of the field perturbations over many expansion times 
on super-Hubblc scales. 
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C. The adiabatic and isocurvature perturbations 

To follow the coupled evolution of the perturbations on large scales (k -C aH) it is more useful to rotate the field 
basis into adiabatic and entropy perturbations, 8a and 8s. Following Rcf. [9], the local rotation is given by 



S Z)=^(lt) where S = ( ™° ~ ^ ) (19) 



8s j V 8<t>2 j \ sin 9 cost 

and the local rotation angle 9, which is now a function of time, is given by 

tan6> = ^. (20) 

The adiabatic and isocurvature perturbations arc then related to the decoupled perturbations Uj by 



aSa \ _ ( cos(9 - 9) - sin(6 -6) \ ( v x 
aSs ) ~ I sin(6 - 6) cos(9 -9) J I v 2 



(21) 



The eigenvalues of U, given in Eq. (11) can be written more compactly in terms of the a and s slow-roll parameters 
(see appendix A) 

Ai,a * \ 



4e - (r ha + r hs ) ± yjw 2 + 4rf as , (22) 

^ L J * 

where 

w = 2e - (r? CTer - T] ss ) , (23) 
and the global rotation angle defined in Eq. (9) simplifies to 

w sin 29 — 2r\ (ys cos 29 



tan 29 = 



w cos 29 + 2r\ as sin 29 



(24) 



From the above equation it is easy to see that the two rotations, by and 9, are the same if and only if = 0, 
i.e. in the special case of a straight background trajectory since [3] 9 ~ —Hrj as . 
Since v% and v 2 are uncorrelated we have 

« 2 (IH 2 > = \((v 2 1 ) + (v 2 2 )) + \cos2(e-9)((vj)-(v 2 2 )) , (25) 
a 2 (SaSs*) = i sin 2(9 - 0)((w 2 ) - («!)), (26) 

« 2 (l^| 2 > = ^((« 2 > + fe 2 »-^cos2(e-0)(( U 2 )-( U2 2 )), (27) 
and from (20) and (24) the trigonometric terms can be written in terms of the slow-roll parameters 

cos2(9-#) = ^ = + 0(e,^,^ s ,?7 sa ), sin2(e-#) = > - = + 0{e, Vaa , r) as ,r)ss) ■ (28) 

V™ 2 + V™ 2 + 

A key observation is that 8a and 8s are correlated at Hubble exit at first order in slow roll if the background trajectory 
is curved. Using the above relationships and Eq. (18) we can calculate the power spectra and cross-correlation at 
Hubble exit 



V„(k) ~ I^J (\ + {-2 + QC)e-2Ci laa ) , (29) 

C as *(k) ~ -2C^ S f , (30) 

P s *(fc) ~ (1 + (-2 + 2C)e - 2Cr] ss ) , (31) 

where here, and in the following, the slow-roll parameters are to be evaluated at Hubble exit, and we define the 
autocorrelation V x = C xx and the cross-correlation 

Att k 3 

C xy S 3 (k-k') ee ^(^•(k')) . (32) 



5 



III. OBSERVABLES AND FINAL POWER SPECTRA 

The comoving curvature perturbation during inflation is given in terms of the adiabatic field perturbations in the 
spatially flat gauge, by [9] 

TT 

K = —8a. (33) 
a 

Similarly a dimcnsionless isocurvature perturbation during inflation can be defined as [2] 

S = ^-5s . (34) 
a 

A convenient dimcnsionless measure of the correlation angle A is 

cosA = ^7^171- ( 35 ) 

It is straightforward to convert the power spectra and cross-correlation at Hubble exit, given in Eqs. (29-31), to 
the dimcnsionless variables 

Vn * - (vL) ( i +(- 2 + 6c ) £ - 2c ^)=^ o) ( i + a i)' ( 3e ) 



Cks* - ~2C Vas = rra 2 , (37) 



V s * - (J^) (l + (-2 + 2C)e-2Cr ]ss )^vi 0) (l + a 3 ), 



(38) 



cos A* ~ -2Cri as , (39) 



where we have defined 

Vi 0) = ( . 

^(T*27T J \ 2lT 

and the first-order slow-roll corrections are given by 

oi = -2e + 2C(3e- Vaa ) 7 a 2 = -2Crj as , a 3 = -2c + 2C{e - Vss ) . (41) 

A. Super-Hubble evolution 

In order to calculate the primordial curvature and isocurvature perturbations some time after inflation has ended, 
which can be constrained by observations today we need to model the evolution of 1Z and S on large scales by 
introducing the transfer functions 



K \ f 1 T ns \ n 
S I T S s )\S 



(42) 



T-rtz = 1 because in the absence of isocurvature modes the adiabatic perturbation is conserved and Tgu = because 
the adiabatic perturbation cannot act as a source to the isocurvature perturbation. 

On the other hand, Tn$ and T$s , which paramcterise the effect of entropy perturbations during inflation upon the 
primordial curvature and isocurvature perturbations, depend upon the full evolution on supcr-Hubble scales, both 
during inflation, and afterwards. Quite generally, we have 

K = a{t)HS, S = P{t)HS , (43) 

and thus we can write [3] 

T ss {U,t) = cxp (3(t')H(t')dt'^j , T ns = J a{t')T ss (U,t')H{t')dt'. (44) 
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The transfer functions contain an implicit scale dependence, through the dependence of the Hubble-exit time, £*, 
which varies with scale, k = (aH)*. The time dependence, calculated using the Leibniz integral rule, is 

1 dT S s a rp 1 dT ns 

H*-du = - /3 * Tss > H,^7 = - a *-^ Tns - (45) 

To calculate the spectral indices of the primordial power spectra to second order in slow roll it is thus necessary to 
also calculate a* and j3* to second order. 

On large scales the time derivative of the curvature perturbation is [19] 

K = 29S. (46) 

Differentiating Eq. (20) for tan 8 gives 

= if sin0cos0(<5f -8"), (47) 

where 8f are slow-roll parameters defined in appendix A. Then using appendix A to write 8f in terms of the potential 
slow-roll parameters we find 

8 2 2 \ 2 

"2 + -£ - -T)<r<7 - ^Vss j ifrs - ~Cas ■ ( 48 ) 

As can be easily seen from (46), cy* = in the case of a straight background trajectory in field space, and the curvature 
perturbation in that case is constant during inflation on supcr-Hubblc scales. 
To calculate /?* we start with the evolution equation of 8s [9] 

8s + 3H8s + (^ + V„+^)s.= i^^ (49) 

where the Bardeen potential \& ~ 4irG(& / H)8a during slow-roll on large scales [9]. Therefore we can drop the term 
on the right hand side of (49) as well as the (k 2 /a 2 )8s term on large scales. Rewriting the equation in terms of S we 
then find 



2^ 



^+{ 3H + 7H)S+[Vss + 39^ + ^H^ + -- w - l ^- w j jS = 0, (50) 
so, to lowest order in slow roll, we have 

S-H(-2e + Vaa - Vss )S. (51) 

We are able to calculate the next order corrections by taking the derivative of this and substituting it back into (50). 
After some calculation, we find 

0* = I (-6e + 3^ CT - 3r] S8 + 126 2 + ?f aa - rf ss - 10er laa + 2er /ss + £ aa - g ss ) . (52) 

B. Final Power Spectra and spectral indices 

The primordial curvature perturbation, during the radiation-dominated era some time after inflation has ended, is 
given on large scales by 

TZ = ^ + ^. (53) 
P 

It is this curvature perturbation that, for example, produces large-scale anisotropies in the cosmic microwave back- 
ground. The conventional definition of the primordial isocurvature matter perturbation is given relative to the 
radiation density by 

s = H (Sp I n_8_p 1 \ 
V Prn Pi ) 
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Prom the definition of the power spectra (32) and the transfer functions (42) it follows that 

V n ~ Vn* + 2TnsCns*+T^ s V s *, (55) 
Cns — T ss Cks* + TssT-jzs'Ps* > (56) 
V S ~ Ti s Vs*, (57) 

and substituting in the power spectra at Hubble exit (36-38) we see that the primordial power spectra are 

V n - Pi° ) {l + Tl s + a 1 + 2T KS a 2 +Tl s a 3 ), (58) 
Cns - V\ ,T S s{Tns + a>2+TKsa 3 ), (59) 
V s ~ ri 0) T§ s (l + a 3 ), (60) 

where oi, &2 and 03 are the corrections from 7-* for the initial power spectra as defined in Eqs. (40) and (41). 

Note that the next-to-leading order corrections to the power spectra (58-60) can be related to the tilts by the 
relation given in Ref. [4] 

Vx-Vf {l + n T -Cn x ), (61) 

where the subscript X means that that the relation is true for the adiabatic, isocurvature and tensor power spectra 
and the cross correlation, and Vx 1S the power spectra at lowest order. The tilt of the gravitational wave power 
spectrum, tit, is defined in section HID 

The scale dependence depends on both the initial power spectra and the transfer functions, we can replace the 
dependence on the transfer functions with the observable correlation angle cos A, which at lowest order satisfies 



;osA (») = -i. (62) 



For details of the calculation and definitions of cos A, sin A, tan A see appendix C. 

The final scalar tilts, defined as nx = dhxVx/dhxk, up to second order in slow-roll parameters are 1 

n n = -(6 - 4cos 2 A)e + 2sm 2 Ar) rTcr + 4 sin A cos A?? CTS + 2cos 2 A7? ss (63) 

10 \ 2 /2 \ 

-— -4cos 2 A +C(24- 16cos 2 A) e 2 + -sin 2 A?^ CT + - cos 2 A +4C(1 - 2cos 2 A) r/ 2 s 



cos 2 A 77 2 s + ( -2 + — cos 2 A + C (-16 + 12 cos 2 A ) J ei] aa - - sin A cos A (1 + 6C)er] as 

2 4 4 

— - cos A (—1 + 6C)erj ss + — sin A cos A (1 — 3C)r/ crcr ?7 (TS + - sin A cos A (1 + ZC)r} as r} ss 

+ 1 sin 2 A (1 + 3C)£ aa + I sin A cos A (1 + 3C) £a S + \ cos 2 A (1 + 3C) £ ss , 
n c = -2e + 2 tan A i) as + 2r) ss (64) 
~ + SCj e 2 + I (1 - 6C) 77 2 s + ^ql + I (l - «W - |(1 + 6C) tan Ae^ + H(l - 6C)er? sa 

+ |(1 - 3C) tanAr?^^ + |(1 + 3C) tan A W7„ + |(1 + 3C) tan A £ as + |(1 + 3C)£ 2 SS , 

n 5 Si -2e + 2?? sa (65) 

22 \ 2 2 8 / S \ 2 2 

-y + 8Cj e 2 + - (1 - 6C) r, 2 s + - ?/s 2 s + - f 1 - -Cj a ha + - (1 - 6C) e %s + - (1 + 3C) £ 2 SS ■ 

The running of the spectral index is defined as ax = dn x /dlnk. Assuming that the power law approximation for 
Vu,Cks and Vs are valid, the running will be second order in the slow roll approximation because the time derivatives 



1 In this notation a scale-invariant (Harrison-Zcl'dovich) spectrum corresponds to nji = 0. 
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(or cquivclantly In k derivatives) of the first-order slow-roll parameters are second order. To leading order in slow roll 
(see also [4]), 

a n ^ 8 (-3 + 4 cos 2 A - 2 cos 4 A )e 2 + 4 sin 2 A cos 2 A rfc a + 4(1 - 4 sin 2 A cos 2 A )j] 2 as + 4 sin 2 A cos 2 A n 2 s (66) 
+4(4 - 7 cos 2 A + 4 cos 4 A )^ aa + 32 sin 3 A cos A en as + 4 cos 2 A (5-4 cos 2 A )e?y ss 
—8 sin A cos A (1 — 2 cos 2 A )ri aa rj as — 8 sin 2 A cos 2 A T) aa r\ ss + 8 sin A cos A (1 — 2 cos 2 A )ij as rj ss 
-2 sin 2 A £ 2 CTCT - 4 sin A cos A £ 2 CTS - 2 cos 2 A £ ss , 
a c S -8e 2 + 4(1 - tan 2 A > ? 2 S + 4e?j CTCT + 4er? sa + 4 tan A n ar7 r) as - 4 tan A n aa ij aa - 2 tan A £ 2 CTS - 2£ M . ( 67 ) 
a 5 £* -8e 2 + 4^. + 4er W + 4e? ?ss - 2£ 2 SS . (68) 

It would be a straightforward but long calculation to include the next to leading terms also, since all that is required 
is to differentiate the tilts including the leading corrections as given above. 

The spectral indices of Vn and V s are both slow roll suppressed while generically ct-ji/s ~ so tne P ower 

spectra are both weakly scale dependent and well approximated by power laws. However the cross-correlation can be 
strongly scale-dependent if tan A is large, specifically if 

tanA?? CTS ~ 0(1). (69) 

In this case the running ac ~ C((tan A rj as ) 2 ) is also large and so C-jis is not well approximated by a power law. 
The problem is that T-jis does not have a power law shape whenever |a*| > \/3*Tjis\ as shown by (45) and then 
Cizs w nl not be close to power law either (56). We can see in this case Cns has an approximate log-dependence on 
wavenumber from (45) and (59). Fortunately this problem only occurs when cos A w 0, i.e. the perturbations are 
nearly uncorrelated and then nc will be nearly unconstrained by observations. Because T-jis has to be small when it 
is not a power law (recalling that tan A ~ 1/T-rs), 'Ptz remains a power law in spite of its leading order dependence 
on T^ s (55). We can see this more explicitly by parameterizing V-r as the sum of two power laws which we discuss 
next. 



C. Alternative parametrisation of the power spectra 



A different way to write the primordial power spectra is to split the adiabatic power spectrum into a part generated 
by the inflationary adiabatic perturbations and a second adiabatic perturbation generated from the inflationary 
entropy perturbation [20], (see also [2, 21]), 

Vn = Al[-) +Al[-\ , (70) 



k J s \k 

k 

B 2 [f) n "\ (72) 



where kg is the pivot scale. The amplitude of the primordial curvature perturbation spectra are given in terms of the 
power spectra at Hubble exit and the primordial transfer functions (see Eq. (55)) by 

A 2 r = [Vn*) ko , Al = [T ns V s * + 2T nS Cns*] ko ■ (73) 

A 2 and A 2 can also be simply written in terms of the correlation angle, (35), 

A 2 = [Pn sin 2 A ] kQ , A 2 S = [P n cos 2 A ] kQ . (74) 

Only 3 of the four tilts are independent because n cm = (n ac i2 + n; so )/2. C-ns and Vs are the same as in the standard 
notation so n COI = nc , rii so = rig and n, ac i2 = 2nc — rig . n a di is the tilt of the adiabatic perturbations at the Hubble-exit 



9 



time of ko ■ The four spectral indices are therefore 

f 10 \ 2 2 /2 \ 

nadi = -6e + 2r? CTCT + ( — + 24Cj e 2 + -r^ + -(1 + 6C)rf as - (2 + 16C)e^ CT + f - + 2C\ & atr , (75) 

(22 \ 2 2 /8 \ 

— + 8Cj e 2 + -(1 - 6C)?7 2 S + -r /s 2 s + ( - - 4Cj e ?W (76) 

4 \ /2 \ 4 4 

- +8CJ tanAe^ + I - -AC J er? ss + -(1 - 3C) tan A i) aa j\ aa + -(1 + 3C) tan Ar) as ri ss 

+|(1 + 3C) tan A£L S + |(1 + 3C)£„ , 

while n cor = nc and rii so = rig are written explicitly in (64) and (65) respectively. 

In the case when tan A rj as is large n ac i2 is large and so is a n ad 2 = % a c — a s ■ Hence the second term of the primordial 
adiabatic power spectrum (70) is not well parameterised by power law. However V-r, can still be well approximated 
by a single power law because A 2 <§C A 2 , from (74) and the requirement that | tan A | 3> 1. 



D. Gravitational Waves 



Scalar and tensor perturbations are decoupled, so the gravitational wave power spectrum is the same as in the 
single- field result, and the amplitude of gravitational waves remains frozen- in on large scales after Hubble exit during 
inflation, [5], 



V T = Vt* ^ 64ttG 



2^ 



(l + 2(-l + C)e) . 



The tilt and running at second order in slow roll can easily be calculated from this power spectrum to give 



tit = — 2e 



11 



-AC 



2C J 7/. 



ax 



-Se 2 + Aeri a 



The tensor-scalar ratio 2 at Hubble exit is the same as in the single-field case 



16e 



- + AC e + - + 2C) r), 



(77) 

(78) 
(79) 

(80) 



However some time after Hubble exit the scalar curvature perturbation may have evolved due to the effect of non- 
adiabatic perturbations. From the definition (35), the final power spectra (58-60) and the curvature power spectrum 
at Hubble exit (36) it follows that 



Vn* = V n sin 2 A 
and thus in the radiation-dominated era we have 



We sin 2 A 



1 



- + AC e + - + 2C U 



(81) 



(82) 



Note that the only difference from the single field tensor-scalar relation is the addition of the sin 2 A factor, this explains 
why the observational upper bound on r does not provide a direct upper bound on e in the multiple field case. 



2 Note that there are various definitions for the tensor-scalar ratio, we use the definition which satisfies r ~ 16e to first order at Hubble 
exit as used in [21], for example. 
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IV. MODEL DEPENDENT RELATIONS 
A. Straight background trajectory 

In the case of a straight background trajectory in field space (during the time of Hubble exit for observable 
modes) the calculation simplifies considerably since then 9 = 0. From (48) this requires that the slow-roll parameters 
7 hs = Cctcts = 0- I 11 tne case of a symmetric potential [11] the background trajectory will be straight since any 
orthogonal velocities will decay quickly. A straight background trajectory is also the attractor solution of some 
assisted inflation models [13-15]. The inflaton and curvaton scenarios are two classes of inflation models with a 
straight background trajectory which we discuss in the following two subsections. 

There is a consistency relation independent of the gravitational wave background which holds whenever the back- 
ground trajectory is straight and which is valid at all orders in slow roll. From (49) the adiabatic and entropy 
perturbations decouple exactly when 9 = 0. Therefore the rotations performed before and after Hubble exit are 
the same, 6 = 9, as was already shown to lowest order in slow roll from (24) and the adiabatic and isocurvature 
perturbations are uncorrelated at Hubble exit. Equations (56) and (57) with C-jis* — imply 

TssC-rs = T-jzsVs- (83) 

From (46) it follows that a* = and therefore from (45), Tss and T-rs have exactly the same k dependence. Taking 
the derivatives of (83) with respect to wavenumber shows that 

n c = n s , (84) 
a c =a s , (85) 

and these consistency relations are true to all orders in slow roll. We can check explicitly that it holds up to second 
order in slow roll from Eqs. (64) and (65). 

The tilts (63) and (65) simplify somewhat to give 

nn = - (6 - 4 cos 2 A )e + 2 sin 2 A j)„„ + 2 cos 2 A n ss + f- y - 4 cos 2 A + C (24 - 16 cos 2 A ) \ e 2 (86) 

2 2 / 14 \ 2 

+ - sin 2 A nl a + - cos 2 A rf ss + ( -2 + — cos 2 A + C (-16 + 12 cos 2 A ) J tr\ aa - 5 cos 2 A (-1 + 6C*)e?7 ss 

+|(1 + 3C) sin 2 A^ CTCT + |(1 + 3C)cos 2 A£ 2 ss , 
n s = -2e + 2 Vsa + (~f + &j e 2 + ^f ss + Q - 4Cj e Vaa + Q - AC^j e Vss + |(1 + 3C) cos 2 Ae 2 ss , (87) 
and the running is given by 

an = 8(-3 + 4 cos 2 A - 2 cos 4 A )e 2 + 4 sin 2 A cos 2 A rf aa + 4 sin 2 A cos 2 A r? ss + 4(4 - 7 cos 2 A + 4 cos 4 A )tr\ aa (88) 

+4 cos 2 A (5-4 cos 2 A )tr\ ss - 8 sin 2 A cos 2 A r] aa r] ss - 2 sin 2 A ^ aaa - 2 cos 2 A ^ 2 SS , 
a s S -8e 2 + 4ei] aa + 4er) ss - 2^ 2 SS . (89) 

B. Inflaton scenario 

If the background trajectory is straight during inflation then this direction in field space, a, can be identified as 
the inflaton field. Other fields orthogonal to the inflaton are time-independent along the background trajectory, but 
quantum fluctuations in light fields (ri ss < 1) do generate isocurvature perturbations during inflation. If we further 
assume that these isocurvature perturbations remain decoupled from the inflaton and radiation density during and 
after inflation then the curvature perturbations today are purely due to inflaton field perturbations, i.e., T-rs = 0. 
We refer to this as the inflaton scenario. Hence any isocurvature perturbations that survive into the radiation era are 
uncorrelated with the curvature perturbations, Cus = and cos A = 0. 

The tilts in the inflaton scenario reduce to 

n n = -6e + 2r? CTCT + (~ + 24c) e 2 + ^ 2 CT - (2 + 16C) &\ aa + f | + 2Cj & aa , (90) 
ns = -2e + 2n ss + (~ + 8C\ e 2 + 2 -r, 2 ss + (| - 4c\ e Vaa + Q - ic\ e Vss + H(l + 3C)^ 2 S s , (91) 
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and the running is given by 



a n = -24e 2 + 16er) aa - 2& aa , 
a s = -8e 2 + 4e?7 CTCT + 4eq as - 2& 



(92) 
(93) 



The results for the curvature perturbation spectrum in the inflaton case is identical to the standard result for single- 
field inflation [5]. Since the curvature perturbations are frozen in from Hubble exit in this case, and the gravitational 
waves also are frozen in, the tensor-scalar ratio is also unchanged from the single field case 



16e 



1 



o + 2C W 



and the standard single field consistency relations apply, see [22] 



-8riT 



l--n T 



n-R. 



(94) 



(95) 



Note that by differentiating this expression we can obtain an infinite hierarchy of consistency relations at higher order 
in the slow-roll parameters [22], the first of which is [12, 23] 



(XT 



(96) 



C. Curvaton scenario 



In the curvaton scenario [16] the primordial curvature perturbation during the radiation dominated era is generated 
from isocurvature field fluctuations in a curvaton field during inflation. These curvaton perturbations lead to an 
inhomogeneous radiation density after inflation if the curvaton is a weakly coupled, massive scalar field whose energy 
density is non-negligible when it finally decays into radiation some time after inflation. 

If the effect of curvaton perturbations is assumed to dominate over the contribution from the inflaton, our chosen 
normalisation of the dimensionless isocurvature perturbations during inflation (34), implies that |T7js| 3> 1. Thus we 
have cos A = ±1 in Eq. (35). The scalar tilt (63) and running (66) are therefore 

n n = -2e + 2 Vss + (~ + 8C^j e 2 + 1$, + Q - 4Cj e VrTrT + Q - 4Cj e Vss + |(l + 3C)£ SS , (97) 
a n S -8e 2 + Ati laa + Ae Vss - 2&„ . (98) 

If the inflaton perturbation, i.e. the scalar perturbation at Hubble exit, is negligible, then from Eq. (80) so is the 
gravitational wave background. We see that for sin A = Eq. (82) implies r = 0. 

Isocurvature perturbations offer an alternative consistency relation in the curvaton scenario. Not only the pri- 
mordial curvature perturbation but also any residual isocurvature perturbation [24] is due to the curvaton field 
fluctuations during inflation. Thus the primordial curvature and isocurvature perturbations are 100% correlated (or 
anti-correlated), and we must have 

nn = n c — n s , (99) 
an = ac = as ■ (100) 



D. Curvature perturbations from broken symmetries 

In [25] it was proposed that the curvature perturbation seen today may be due to the isocurvature perturbations 
during inflation in a nearly symmetric potential. The isocurvature perturbations are then converted into curvature 
perturbations during instant preheating at the end of inflation [26]. The key parameter of the preheating is the 
minimum distance to the minimum of the potential reached along the inflaton trajectory and perturbations of this 
minimum distance are due to the isocurvature perturbations. 

The nearly symmetric 2-ficld potential considered in [25] is 

nfc,fc) = x(^ + iTi) . ( 101 ) 
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where x is the symmetry breaking parameter which is assumed to satisfy < x <C 1. We will assume the symmetry 
is only weakly broken and so x is of the same order as the slow roll parameters or less. Hence we will neglect terms 
of order xe 2 and x 2 e. 

The values of the slow-roll parameters depend on the angle of the background trajectory which is a free parameter 
of the theory The distance from the minimum of the potential when observable scal es today were leaving the horizon 
during inflation (i.e. about 60 efoldings before the end of inflation) is fixed at a = \f <j>\ + (f>?, — 3Mp;. The <fii and 02 
field values can be related to the differential rotation angle (19), defined by (20), by 

0i ~ crcos6»(l -aisin 2 6>), 2 ~ crsin6>(l + a; cos 2 9) . (102) 

These can be derived by applying slow-roll approximations to the equations of motion 

2 

01 + 3F0i + m 2 x = , 4> 2 + 3ff0 2 + 02 = 0. (103) 

1 + X 

Note that in the case of a symmetric potential, (x = 0), the background trajectory is straight with a polar angle 
coinciding with 9 [11]. 



The potential can be written as 



The slow-roll parameters are therefore 
1 1 



V~ -^-a 2 (l-xsm 2 6). (104) 



4nG a 2 



V™ = e, 7]as ^ -xe sin6 cos 6, r) ss ~ e(l - a; cos 20) , £ CT(r(r = ^% as = £ ass =0. (105) 



If we suppose that the isocurvature perturbation during inflation dominates the primordial curvature perturbation 
we have cos A = 1 and the spectral tilt is 

n n = -2xe cos 29 - — e 2 . (106) 
3 

Note that it is quite possible for the second-order in slow-roll contribution to dominate over the first-order result. 
In the limiting case of a symmetric potential, V = to 2 ct 2 /2, there is no contribution at first order and the tilt is 
nn = — 10e 2 /3 as calculated in [11] using another method. 



V. CONCLUSIONS 



We have calculated the primordial curvature and isocurvature perturbations for an arbitrary model of two field 
inflation including first-order corrections in the slow-roll expansion. We have calculated the power spectra and the 
cross-correlation, including the first order corrections [Eqs. (58), (59) and (60)] which allows us to calculate the 
spectral tilts to second order in slow roll [Eqs. (63), (64) and (65)]. We find as expected that the scale dependence of 
the curvature and isocurvature power spectra are small, being first order in slow roll. However the scale-dependence 
of the cross-correlation is not necessarily small. It becomes large for r\ as > cos A , where the slow-roll parameter 
r] as , defined in appendix A, determines the curvature of the trajectory in field space. Similarly we find that the 
running of the curvature and isocurvature tilt is very small, second order in slow roll, and thus the scale-dependence 
is well-fit by a power-law. However the scale dependence of the cross-correlation is not well-described by a power-law 
if rj as > cos A . 

We use the approach of Gordon et al. [9] to evolve the instantaneous adiabatic and entropy field perturbations on 
super-Hubblc scales. This approach was generalized to an arbitrary number of fields in curved field space by Nibbclink 
and van Tent [4, 10], who also calculated the leading slow- roll corrections to the primordial curvature and isocurvature 
power spectra, making some assumptions about the evolution of the universe after inflation. The isocurvature pertur- 
bations can be evolved independently of the curvature perturbations on large scales, but the curvature perturbation 
can be altered by non-adiabatic perturbations. The primordial power spectra are thus dependent upon two transfer 
functions, T-rs and whose magnitude we leave arbitrary, but whose scale-dependence can be determined in terms 
of the slow-roll parameters at Hubble-exit during inflation [3] . 

However on small scales (at early times during inflation) the coupling between the instantaneous adiabatic and 
entropy modes can become large leaving the initial vacuum state ambiguous in terms of these variables. Instead we 
defined an orthonormal basis in field space in which the interaction between fields is negligible at early times and 
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minimises the interaction for a given mode at Hubble-exit. In this basis the effective mass matrix is diagonaliscd and 
the field perturbations are uncorrelated at Hubble exit at first order in slow roll. But we find that the instantaneous 
adiabatic and entropy field perturbations of Gordon ct al. [9] are correlated at Hubble exit at first order in slow roll, 
cos A * oc rj asi unless the trajectory in field space is a straight line. Correlations between perturbations during 2-field 
inflation was also studied in [18]. 

In [3, 27] it was shown that the well known single-field consistency relation r ~ — 8nr can be generalized to the 
case of two-field inflation to give r — — 8riT sin 2 A . Including the next-order correction terms the two-field consistency 
relation becomes 

r = —8riT sin 2 A 

At leading order in slow roll with two fields, there are four slow- roll parameters (e, r\ ac! , rj as and f] ss ), and two 
transfer functions (Tns and Tss) as well as the energy scale (or Hubble parameter) at Hubble exit during inflation. 
However there are, in principle, eight observables, corresponding to the primordial curvature and isocurvature power 
spectra, their cross-correlation spectrum, the tensor spectrum and the tilts of all four spectra. Thus there is one 
consistency relation [3]. At second-order we have three more slow- roll parameters (£,„ aa , £ 2 CTS and £,% ss )- Although 
£ 2 SS is in general non-zero it does not appear in the calculation of the power spectra, because whenever we take a time 
derivative (or derivative with respect to scale) we are differentiating along the background trajectory, and therefore 
with respect to the adiabatic field, a and not s. But we also have four more observables at this order, corresponding 
to the running of the four spectra. And thus we have one more consistency relation. Of course, there is no guarantee 
that either the primordial isocurvature perturbations or gravitational waves will be large enough to be detected. 

We can differentiate Eq. (107) with respect to wavenumber, to find a higher-order consistency relation in terms of 
the running of the tensor spectrum, see [4] , 

cxt sin 2 A = tit [tit — Tin + cos 2 A (2nc — n$ — tit)] ■ (108) 

We could include the third-order correction to this consistency relation by including the second-order correction term 
in (107), but the result would be extremely long. Since (107) holds on all scales we could repeatedly differentiate it 
to calculate an infinite hierarchy of consistency relations, as done in the single field case [22]. 

In the case of totally correlated curvature and isocurvature perturbations (i.e. cos A = 0) the first consistency 
relation (107) is trivially satisfied since both r = and sin A = 0. This corresponds for example to the curvaton 
scenario discussed in section IV C. However in this case there is an alternative set of consistency relations, rig — nc = 
n-jz at first order and as = ac — ct-jz at second order and so on to give an alternative infinite hierarchy of relations. 

Finally note that observational constraints on the tensor-scalar ratio in the single field case directly constrain e 
to be small, and then the near scale invariance of the scalar tilt also constrains r\ to be small. Hence higher-order 
slow-roll corrections in single field inflation are constrained to be very small. However in two-field inflation the upper 
bound on the tensor-scalar ratio does not provide a direct constrain on e, and the near scale invariance of the adiabatic 
power spectrum only constrains a certain combination of the slow-roll parameters to be small. Hence the higher-order 
slow- roll corrections in two-field inflation can be significant. 
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APPENDIX A: THE SLOW-ROLL PARAMETERS 

Different papers use different definitions of the slow-roll parameters, some of which are equivalent at first order, 
but none are equivalent at second order. The only ones which we shall use are 



, 1 1 

2 sm A 



n-ii 



„ cos~A cos A 
2 — — nc + . 9a n s 



sin^A 



sin 2 A 



(107) 
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First order Hubble 



F H ------- e H +f H fAll 

e - H 2 - 2 IP ~ 22 ' [ ' 

H _ 8ttG 



c " ~ 2 h 2 



First order potential 



1 (V, 



16-kG \ V 



Second order potential 



(A2) 



S? = -4- & . (A4) 



(A5) 



Vas - SirG V ' n °° ~ SttG V ' V " ~ 8ttG V ' ( ' 

^ ~ (8ttG) 2 y 2 ' ~ (87rG) 2 y 2 ' 1 ' 

f 2 VcrssVg si 1 VsssVg ( AR) 

?<TSS ~ (8ttG) 2 V 2 ' ?sss ~ (8ttG) 2 V 2 [ 1 

Note that the superscript 2 refers to the quantity being a second order slow-roll parameter. 

Slow-roll parameters defined in terms of the Hubble parameter can be related to the potential slow-roll parameters 

by 



where 



e H = e--e 2 + -a ha , (A9) 

€ = -C + + ^ ~ 3«7«r + + \nls + , (A10) 

Sf = -e H + A I + Uj(Aj - e) - ieAj - , (All) 



sin cos 6* 

Ai = 7711 H -7712 , A 2 = 7722 + ~. 77^712 , ( A12 ) 

cos 6< sin # 

„ , *ha n cos 2 0-sin 2 , Aio , 

Al + A 2 = Vn + 2) + ? /22 = 2r)a* + Vas — D ^— ■ (A13) 

sm cos # sm 6* cos a 



APPENDIX B: DERIVATIVES 



Note that all slow-roll parameters can only be taken to be constant at the order at which they are defined, i.e., the 
time derivative of a first-order slow roll parameter is second order. 

^£ = 26(26-^), (Bl) 

-^77™ = 2ei laa - 2r, 2 as - £ arT , (B2) 

-^■i] as = 2ei]as + Vas(V<™ - rise) - Cas , ( B 3) 

^7) ss = 2e Vss + 2?7 2 s - C ss , (B4) 

l^ = -4e 2 -2r7 2 s +4e77 ffff -£ 2 CTCT - (B5) 
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All of the relations above will also be true to the same order if jjj^ is replaced by d ^ nk at k = aH. 

The derivative of e H is also required to third order for calculating the spectral tilt of vi^ (40) to second order in 
slow roll, and for calculating /3* (52) to second order in slow roll, 

d\ne H /„ 7 „ 2 1 2 1 9 1 2 \ 



d\nk V ' 3 3 3 ' CTS 3 



APPENDIX C: CALCULATING THE TILT 



In order to calculate the spectral indices of the power spectra (58,59,60) we need to calculate the scale dependence 

of Pi°\ the transfer functions and a\,a 2 , 0,3, (41), all up to second order in slow roll. 
Wc define 

n ^ = "TbgT" (C1) 

where is defined by (40) and split the result into first and second order parts 

n (oj = ~ 6e + 2r l<y<yi n [o) = y fi2 + \~nla + 1 7 lls - 6e?7a<7 + ^Caa ■ (C2) 

The scale dependence of the transfer functions (45), where the functions a* and j3* are given by (48,52) respectively 
and using dlogfc ~ H(l — e H )dt* are 



d\ogT ss 
dlogk 



(i) 



2e - r 1aa + r lss , (C3) 



3lQgTgg (2) „ 2 1 2,1 2 ,7 1 1 2 1 2 

-gj^- = -2e - 3^ + + ^ + - + 3^ , (C4) 

^logT^W 



d log 



= 2e - rjo-CT + 2 tan A 77^ + ?7 SS , (C5) 



^iF~ 2 = ~ 2 ' 2 " 3^ " 4Ctan2Aj ?- + ^ + - tan A (| + 4C ) + ^ ( C6 ) 

+ tan A Q + 2Cj 77^77^ + tan A (| - 2C) Mss - ^g aa + | tan A + • 

Note that the a/, (41) are first order in slow roll, so their derivatives are second order in slow roll and can be easily 
calculated from (B1,B2,B3,B4). 

In order to write the spectral indices in terms of observables we need to be able to relate the transfer functions and 
the correlation angle (35). From (58,59,60) it follows that at zeroth order 

C0S A(°) = — sinA^ = ^^=, tanA(°) = - , (C7) 

and at first order 



cosA (1 > =cosA<°) 
sinA (1) ee sinA ( °) 



1 . 2a surA 1 . o. 
— a\ sin A + ci2 7 — I — 03 sm A 

2 cos A 2 

\-a\ cos 2 A — a 2 sin A cos A — ia 3 cos 2 A 



(C8) 
(C9) 



Another useful equation for Tug is 

1 



= tanA (0) - tan A (1 + C (— 2e + r\ aa - 2rj as tan A - rj ss )) . (CIO) 
ins 
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To find the spectral tilts we take the log derivatives of (58,59,60) and substitute away the dependence on the transfer 
functions by replacing them with the correlation angle (C7,C8,C9,C10) to find 



n « = +2sinAcosA-=£S., (Cll) 
w alogfc 

- - (2) ^2sinAcosA^St-2(cosAsinA( 1 )+sinAcosA( 1 ))^i (C12) 
a log k d log k 

dT (1) dT {1) 



K — "-(0) 



-2ai sin 3 A cos A — + 2a 2 sin 2 A (sin 2 A - cos 2 A ) — 

d log k d log k 

q . . drX 1 ^ . ?a <^ a i . . da,2 9 A da3 

+2a 3 sin d A cos A - K ^ + snrA - - + 2 sin A cos A ,, — + cos 2 A ■ 



dlogfc dlogfc dlogA; dlogfc 



(1) 

c 


- n (o) 


d\ogk 


diogr^ 

dlog k 




(2) 

- n (0) ~ 


dlog fc 


d\ogT%l 
dlog k 


(1) 

s 


- n (o) 


d log k 




f 


(2) 

- n (0) ~ 


JlogT^ 
dlogfc 


da 3 
dlogk 



(C13) 

a . dlogT^l da 3 ,_ H 

tan A— ^-- a2 tanA 6 + — ^- , C14 
a log k a log /e a log k 

(C15) 
(C16) 

So substituting (C2,C3,C4,C5,C6) in the above six equations we find the tilts as displayed in (63,64,65). 
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